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Abstract. The high latitude translucent molecular cloud L1642 has been mapped in the J=l-0 and J=2-l transitions of 12 CO, 
13 CO and C 18 using the SEST radio telescope. We have analysed the morphology and velocity structure of the cloud using 
the Positive Matrix Factorization (PMF) method. The results show that LI 642 is composed of a main structure at radial velocity 
0.2kms _1 while the higher velocity components at ^0.5 and 1.0 km s -1 form an incomplete ring around it, suggesting an 
expanding shell structure. Fainter emission extends to the north with a still higher velocity of up to 1.6 km s -1 . Such a velocity 
structure suggests an elongated morphology in the line of sight direction. The physical properties of the cloud have been 
investigated assuming LTE conditions, but non-LTE radiative transfer models are also constructed for the 13 CO observations. 
We confirm that LI 642 follows an r _1 density distribution in its outer parts while the distribution is considerably flatter in the 
core. The cloud is close to virial equilibrium. 

In an Appendix the PMF results are compared with the view obtained through the analysis of channel maps and by the use of 
Principal Component Analysis (PC A). Both PMF and PC A present the observations as a linear combination of basic spectral 
shapes that are extracted from the data. Comparison of the methods shows that the PMF method in particular is able to produce 
a presentation of the complex velocity that is both compact and easily interpreted. 
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^ ■ 1. Introduction wavelength data (Liljestrom and Mattila 1988; Taylor et al. 

I t t , . , , r [19821 Laureiis et al.fl%7l Liliestrom fTWTl L1642 appears as a 

High latitude translucent molecular clouds are ideal targets for T , . T , T . . . . , . U1 UT 

, t r t . . . rr r - cool and quiescent cloud. It is associated with a much larger HI 

the study of the interactions between different phases of the , A , AO , .. u , . .. 

r. t . t . 1 t . t t r t . cloud (over 4 ) with cometary structure. This is also seen m the 

interstellar medium which are subiected to the general inter- Tri AC , nn T , + ■, ■ A . T . , u ^ , 

J fc IRAS 100 /im image. The tail is perpendicular to the Galactic 

stellar radiation field. The visual extinction of these clouds is T A . , . A , u , 

plane, and points towards the plane, 

modest (Ay ~ 1-3 mag) which places them between dark and ^, TnAC . mAon/noci/im a 

i i * -t . . rr ,i i lt i The IRAS point sources IRAS 04325-1419 and 

dilluse clouds. As they are minimally affected by embedded T n\c n^ni/nn i T u.n ^ +u i + + 

J t . . /i i . i i . i IRAS 04237-1419 make L1642 one of the closest star 

star-formation the external radiation field is the dominant heat- € . , A m U T . . u . , u 

t t . t t .1.1 r i forming clouds. The source colours suggest that they are 

mg source. Furthermore, due to the high galactic latitude of the ^ ^ . . T . , rrt A n A ~~ n A * n i TT u- 

t & t . ' t . ii r i T-Taun stars. In particular, IRAS 04237-1419 shows Herbig- 

clouds, observations are seldom contaminated by foreground or TT u _ . /0 A „ ^ 1 i ino ^ i . .* <-> t t . i . 

t ™x . Haro characteristics (Sandell et al. 1987 Reipurth & Heatncote 

background emission. The cloud LI 642 (=MBM20) is a proto- inrkrk i TT ^ ^ i /r \ 

ri . i r i i T . i i i . t . 1990). However, the low luminosity of the source (L < 0.6 L©) 

type of this class of clouds. It is located at galactic coordinates ±u < * ™ u i ^ i • a 

7 7 ^ mi i • • i TT i means that it will have only a very local influence. 

/ =2 10. 9 , b=- 36.5 . The distance determination by Hearty etal. T ^. , . fT 

j— — — i ' 11^^ ^ i v In this paper, we present an analysis of LI 642, making use 

(2000) places it between 112 and 160 pc. Based on X-ray mea- r , ... , 

1 t r lt T i^„^. .1. of new CO observations with a wider map area and/or denser 

surements, Kuntz et al. ( 1997) suggested that LI 642 is withm . , - . x1 . » ^ -. . . . r , 

_ i — 7^ 11111 , . , . , . grid than in the previous studies. After the description of ob- 

or close to the edge of the local bubble, which in this direc- . . rn o . m- j , j , ^ i r 

^ . servations m Sect. |2| Sect. |3Jis devoted to the analysis of mor- 



tion is at about 140 pc according to Sfeir et al. ( 1999). This is 
in good agreement with the direct measurement and we adopt 
140pc for the distance of LI 642. Based on previous multi- 



phology and kinematics using the Positive Matrix Factorization 
method. The analysis is also carried out using another multi- 
variate method, the Principal Component Analysis (Murtagh & 
Send offprint requests to: D. Russeil Heck 1987 ) and by channel maps. A comparison of the analy- 

* Based on observations collected at the European Southern sis methods and some comments on their relative strengths and 
Observatory weaknesses are given in an Appendix. In Sect. |4| we discuss 
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Fig. 1. Averaged spectra observed towards the centre of LI 642 
(04h35ml0.1s, -14°15'58" (J2000)). The solid, dotted and 
dashed histograms are respectively the 12 CO, 13 CO and C 18 
profiles. 

the correlations between the observed lines, and in Sect. |3 the 
physical properties of the cloud are derived. Finally, in Sect. 
we present our conclusions. 

2. Observations and data reduction 

Molecular line observations were made by K.L. in August 
1999. The transitions J=l-0 and J=2-l of 12 CO, 13 CO and 
C 18 were observed with the Swedish-ESO Submillimetre 
Telescope (SEST). The receivers were cooled Schottky diode 
mixers operated in 80-1 16 GHz and 216-245 GHz bands. The 
channel separation of the AOS spectrometer was O.llkms -1 
for the J=l-0 transition and 0.06 km s -1 for J=2-l transi- 
tion (Booth et al. 1989). Observations were made in frequency 
switching mode. The calibration was done using the chopper- 
wheel method (Ulich & Haas 11976ft and checked by observ- 
ing the Orion KL nebula. Antenna temperatures (T£) were 
converted to radiation temperatures (T^) using the Moon effi- 
ciency. The half power beam widths are 45" and 23" for the 1-0 
and 2-1 transitions, respectively. The grid spacing was 3' . The 
grid positions were planned to be identical to the map positions 
of our 200 /im ISOPHOT map (Lehtinen et al., in preparation). 

Fig. [l] shows examples of observed line profiles. Typical 
noise levels were AT rms =0.06 K for the C 18 0(2-l) spectra 
and ~0.15 K for all the other transitions. 
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Fig. 2. Contour maps of the 12 CO line area. The square and 
the triangle indicate the position of IRAS 04327-1419 and 
IRAS 04325-1419, respectively. The two numbers given for 
each map are the isocontour starting values and the steps in 
Kkms -1 . The centre position (0,0) is 04h35m30s, -14°24'28" 
(J2000). The arrow indicates the direction of the cloud tail seen 
in HI and IRAS maps 



3. Morphological and kinematic analysis 

The observed line area maps are presented in Figs. 12 141 We 
utilize the multi-line information to probe the cloud morphol- 
ogy and kinematics from the dense central parts, as traced by 
the C 18 0, to the diffuse outer parts, as traced by 12 CO. The 
conventional channel maps as well as the results from princi- 
pal component analysis (PCA; see e.g. Murtagh & Heck M1987l 
are given in Appendix[A] Here we study the cloud structure by 
using Positive Matrix Factorization (PMF), since comparison 
with the other methods showed that this method gives the most 
compact presentation for the complex cloud structure. The dif- 
ferent analysis methods were not only used to study this partic- 
ular cloud, Lynds 1642, but also to evaluate the relative merits 
of the methods. These are discussed in Appendix lA.3l 
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Fig. 3. As in FigEI but for 13 CO 
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Fig. 4. As in Fig|2| but for C 18 



3. 1. Positive Matrix Factorization 

The Positive Matrix Factorization (PMF) method assumes the 
spectra to be the sum of a few basic velocity components. The 
line decomposition computes the shape and intensity of each 
component imposing positivity to their profile. The application 
of PMF to the analysis of molecular lines was largely discussed 
by Juvela et al.( 1996 ) and we repeat here only some of the main 
points. Let X be the data matrix and E a matrix consisting 
of estimated standard deviations of the elements of X. Both 
matrices have dimensions nxm. Given matrices X and E and 
the rank of factorisation, r, PMF calculates a factorisation 

X = GF + E,G ik >0,F kj >0 (1) 

with i=l,...,n, k=l,...,r and j=l,...,m). E is the residual ma- 
trix and factor matrices G and F are selected to minimise its 
weighted norm 

The rank r equals the number of spectral components used in 
approximating the observed data. 

There are two essential improvements over PCA analysis 
(see Appendix \A.2l . Firstly, the elements of the factor matri- 
ces, G and F are required to be positive. In this case, this is 
a physically meaningful assumption. The observed spectra are 
decomposed into positive spectral components and a direct in- 
terpretation of the results will be easier than in the case of PCA, 
where profiles may contain as many negative as positive values. 
Secondly, the residuals are weighted according to individual 
error estimates of the observations. This makes it possible to 
use all observations without the risk that results will be unduly 
affected by individual, low signal-to-noise ratio spectra. Like 
PCA, the PMF method regards the observations as a sum of 
emission components at fixed radial velocities. In the case of 
continuous velocity gradients, this is not an exact description 
of the data, and this must be kept in mind when interpreting the 
results. 

PMF requires the selection of the number of components, 
r, and a suitable value of r could be determined e.g. based 
on the eigenvalues of matrix X. In practise, it is better to try 
different r- values and select the best value based on the de- 
values obtained and the appearance of the components and the 
residual maps. The PMF method was run independently on the 
13 CO(1-0), 13 CO(2-l), 12 CO(1-0), 12 CO(2-l), C 18 O(l-0) and 
C 18 0(2-1) datasets to extract kinematic information for LI 642. 
Factorisations were computed for ranks r=l to r=5 using the 
Multilinear Engine program (Paatero 1999). Results of factori- 
sations with r = 4 are shown in Figs. The figures show 
the computed spectral profiles of the components for J =1-0 
lines and their spatial distribution over the mapped area. The 
results for the J=2-l transitions were very similar to those of 
J=l-0. By trying decomposition with 1 to 5 components, we 
found that, in 12 CO and 13 CO, 4 components are enough to de- 
scribe the main kinematic structures. For C 18 0, only three sig- 
nificant components were found, but 4 components are shown 
for the sake of completeness. 
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Fig. 5. Result of the PMF factorisation of 12 CO(1-0) observations with four spectral components. The spectral profiles of the 
components are shown in the first frame. The other frames show the intensity distribution of the components over the observed 
map. The component maps are in the order of increasing radial velocity of the component 



3.2. ClOUd kinematics and environmental effects are not pure gaussians, but exhibit small shoulders and slightly 

flattened tops (see Fig. [TJ. The possible causes for such line 
We analyse the kinematic structure of L1642 using the PMF asymmetries are: ejection, contraction, rotation, self absorp- 
results shown in Figs. 15171 The maps of the kinematic compo- 
nents show an intricate structure. Most of observed CO lines 
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Fig. 6. Result of the PMF factorisation of 13 CO(1-0) observations with four components 



tion, saturation, or superposition of several kinematically sepa- 3.2.1. 12 CO lines 
rate cloud fragments along the line of sight. 

With r = 4, the first component of PMF decomposition traces 
the central part of LI 642. The second one indicates strong 
emission south of the centre, extending to the south west 
edge of the map. This is still visible in the third component. 
However, the last two components mostly trace emission from 
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Fig. 7. Result of the PMF factorisation of C 18 O(l-0) observations with four components. The last component is mostly noise 
and is drawn with dashed line 

the northern part, indicating a velocity gradient towards the kms -1 . With more components (r=5), the Q-value improves 

north, or possibly the presence of another cloud component further by more than 15%, and decompositions show more de- 

at a higher radial velocity. The decomposition shows patchy tails of the velocity structure in the northern part. The general 

emission at ~ 1 .25 km s _1 , most notably the prominent clump kinematic structure is clear, but it cannot be interpreted directly 

at (-20 / ,+20 / ) embedded in more extended emission at ~0.4 e.g. as simple rotation, contraction or expansion. The smallest- 
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Fig. 8. Observed 12 CO(1-0) and 13 CO(1-0) spectra (upper 
frames) and the residuals (lower frames) of the PMF fit at po- 
sitions (0', 90 and (-21',21'). The 12 CO spectra are drawn with 
thicker lines 

velocity component 1, centred at ~0.2kms -1 , delineates the 
maximum CO emission core of the cloud. At higher radial ve- 
locities (components 2 and 3) emission in both northern and 
southern parts of the maps can be seen. The highest-velocity 
component 4 is only present in the extreme North. 

Its velocity difference of ~1.3kms _1 with respect to the 
cloud core corresponds to a distance of ~1 pc travelled in 10 6 
years, a typical timescale for cloud dynamics. 

The residual maps of PMF factorizations with r=4 (and 
with r=5) are very similar for the two 12 CO transitions, and 
both show several clear clumps, most notably at positions (0', 
9') and (-21', 2Y). There is also a wider area of large residuals 
between these positions, at around (-10', +100, i-e. around the 
position of the source IRAS 04325-1419, which is associated 
with an outflow (Liljestrom et al. 1989). 

The study of residuals is a very powerful tool for finding 
deviating spectra. By subtracting the decomposition from the 
observation matrix one can not only identify deviating spectra, 
but also see exactly how the profiles differ. For example, two 
distinct structures were observed in the PMF residual maps of 
both 12 CO(1-0) and 12 CO(2-l). Around position (0', 90 the 
spectra show clear blueshifted excess (v ~-0.5kms _1 ), while 
around position (-21', 210 the excess is as clearly redshifted 
(v ~+1.3 km s _1 ) (see FigEJ. These residual features are quite 
strong, with peak 12 CO(1-0) antenna temperature above 1 K. 
Corresponding features can even be detected in channel maps. 
However, the PMF analysis indicates that the emission seen 
in the -0.5 km s -1 channel is not a normal tail of the stronger 
emission at higher velocities. The spectra at the given position 
have distinctly different profiles. 

3.2.2. 13 CO lines 

The 13 CO emission is less extended, but it exhibits similar 
kinematics as 12 CO. Figure shows the 13 CO(1-0) decompo- 
sition with r=4. The maps for the transition J=2-l are again 
very similar. The first component is concentrated in the cloud 
centre, while the last one traces emission in the north, in partic- 
ular the two clumps at (-20', +200 and (-5', +150, already seen 



in the 12 CO maps. On the other hand, the ridge of emission 
south of the centre is weak in 13 CO. For 13 CO(1-0), a change 
from r=l to r=2 decreases Q by 35%, while the drop is only 
15% in the case of 13 CO(2-l). A similar difference is noted for 
higher r- values. This is caused by the difference in the signal- 
to-noise ratio, and with r > 3, some noise ripples appear in 
the 13 CO(2-l) components. As in the case of 12 CO residual 
maps, the 13 CO(1-0) maps show significant residual emission, 
e.g. around position(0 / , +90 (see Fig.[8|). 

3.2.3. C 18 lines 

C 18 probes the deepest and densest parts of LI 642. Compared 
with the previous lines, it is spatially much more concen- 
trated and is restricted to a velocity range between -0.3 
and 0.5 km s -1 . The PMF 1 -component decomposition of the 
C 18 O(l-0) shows that the emission has a mean velocity of 
0.2 km s _1 and a width of 0.47 km s _1 . We show in Fig.0the 
decomposition with r = 4 for completeness, even if decom- 
positions with r > 3 have some irregular components that are 
caused mostly by noise. For the two transitions, the maps of 
the first three components agree with each other and these are 
therefore likely to be real emission features. At low radial ve- 
locities, a core located at (-4', +90 is the main feature. At in- 
termediate velocities, two separate clumps appear at (+3', +50 
and (-4', +60- The third component shows some emission north 
of the core, and, in the maps of both transitions, weak clumps 
are seen at (-9', +120 and (+7', +90- The last components have 
completely different distributions and spectral profiles and, and 
therefore are unlikely to be related to actual cloud structure. 
Like residual maps, these 'noise factors' can be used to check 
for systematic errors or artifacts in data. For example, in the 
case of C 18 0(2-1), the fourth spectral component is almost flat 
and positive, suggesting bad baseline fits in a few spectra. The 
two IRAS sources do not coincide with the main peak of the 
C 18 emission, but are located inside the western clump seen 
in the third component map, i.e. around velocity ~0.3 km s _1 . 

In addition to the clumps already mentioned, there are sep- 
arate patches of emission that are seen in the second and the 
third component maps. They can be either independent clumps 
or parts of a more continuous density distribution containing 
gas motions along the line of sight. The latter hypothesis seems 
more likely, because the features do not have a clear counter- 
part on peak intensity and line area maps. 

3.2.4. Summary of kinematics and environmental 
effects 

To summarise, PMF (with increasing r) provides a hierarchical 
view of the cloud kinematics. Firstly, the main cloud is distin- 
guished from the emissions in the north. Secondly, kinematic 
substructures appear in the northern part, and finally, kinematic 
substructures are identified inside the central area of LI 642. 
The northern part is likely to belong to LI 642, since the ve- 
locity difference is less than the typical velocity dispersion of 
molecular clouds (the typical cloud-to-cloud velocity disper- 
sion is 7±1.5kms _1 , according to Blitz and Fich, 1983|. The 
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Fig. 9. The spatial distributions of the first two PMF compo- 
nents in the case of the 12 CO(1-0) and the 13 CO(1-0) observa- 
tions. The blue and red contours correspond to the first and the 
second PMF component, respectively, as shown in Figs. |3 and 
|6j The contour levels are 40, 50, and 55 for the first and 20, 25, 
and 30 for the second 12 CO component, and 14, 19, and 24 for 
the first and 8.5 and 1 1.5 for the second 13 CO component. The 
units are the same as used in Figs. and |6l 



12 CO observations show a systematic velocity structure with 
the main emission at ~0.2kms _1 (a velocity close to that of 
the C 18 core), and fainter and patchy emission with veloci- 
ties up to 1.6 km s _1 extending to the north and to the south. 

In the PMF decompositions of the 12 CO lines the first com- 
ponent shows the compact core. The second component shows 
almost a ring-like structure at ~0.5 km s _1 higher radial veloc- 
ity. This is illustrated in Fig. [9^ where we have overlay ed the 
contours of the first two PMF components. Similar structure is 
seen in the decomposition of the 13 CO(1-0) observations (see 
Fig-Hb)- This suggests some kind of an expansion. At higher 
velocities, no emission is seen towards the centre of the cloud. 
Therefore, the structure cannot be spherically symmetric, but it 
could be an incomplete expanding shell, the far side of which 
is missing. 



The LI 642 cloud is the head of a large cometary structure 
with a tail extending more than 5 degrees in the north eastern 
direction. The tail is seen very clearly, e.g. in the 100/im IRAS 
map. Furthermore, HI data shows a velocity gradient along the 
tail, with increasing velocity towards the north-east (Taylor et 
al. 119821) . This suggests that the first PMF component repre- 
sents the very head of the whole structure, and gas in the head 
is flowing relative to the rest of the structure, not only towards 
the south-west but also more towards us. The second PMF com- 
ponent would correspond to more diffuse material being left 
behind by the moving cloud core. In the north-eastern corner, 
the spatial separation of the third and fourth PMF components 
indicates a velocity gradient across the tail. This velocity gra- 
dient closely agrees with the HI velocity structure (Lifjestrom 
& Mattila 1988; Taylor et al. fT982b indicating that the atomic 
and molecular gas components are well mixed in the cloud. 

In projection on the sky, LI 642 is located in the direction of 
the edge of the Orion-Eridanus Bubble, and we must consider 
the possibility of their interaction. /From multi-wavelength 
analysis, Reynolds and Odgen \ 19791) have shown that the 
Orion-Eridanus Bubble is an expanding cavity (expansion ve- 
locity ~15 km s _1 ) full of warm ionised gas surrounded by an 
expanding HI shell. X-ray enhancements have been measured 
and attributed to hot gas flowing out from the bubble through a 
chimney or a hole where LI 642 would be located (Heiles et al. 
1999). From HI data (Brown et al. 119951) . the bubble velocities 
range between -40 and +40kms _1 . The largest spatial exten- 
sion occurs between -1 and +8 km s _1 , which is the most prob- 
able range for the systematic velocity of the bubble. LI 642 falls 
in the velocity range of the Orion-Eridanus Bubble, which sug- 
gests that it could have been swept up and hence shaped by the 
bubble. This would favour the hypothesis that the tail is formed 
as a stationary cloud embedded in moving interstellar medium. 
On HI maps (Brown et al. 1995), L1642 appears isolated for 
velocities between -1.5 to 2.5kms _1 . At these velocities, the 
bubble edges delineate an 'empty space' roughly centered on 
LI 642. More precisely, LI 642 is kinematically disconnected 
from the filaments at I ~210° which have velocities less than 
-3kms _1 . 

4. Correlations of the line intensities 

We use the multi-transition information to examine the excita- 
tion conditions of LI 642. There is a good correlation between 
the integrated line intensities of the J =2-1 and J =1-0 tran- 
sitions (Tig.[Tol). The following relations are obtained for the 
line areas: 

W[ 12 CO(2 - 1)] = (0.75 ± 0.06) W[ 12 CO(l - 0)] 

-(0.12 ±0.40) 
W [ 13 CO(2 - 1)] = (0.71 ± 0.01) W[ 13 CO(l - 0)] 

-(0.01 ±0.19) 

(3) 

These relations are in agreement with other high latitude clouds 
(e.g. Ingalls et al. l2Q00b . The 12 CO/ 13 CO ratio is much smaller 
than the typical relative abundance of the species, which is a 
natural consequence of 12 CO being optically thick. 
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Fig. 10. Correlation between 12 CO J =2-1 and J =1-0 line Fig. 12. Excitation temperature T ex (l-0), calculated from 

13 CO/ 12 CO peak antenna temperature ratio, versus radius. 
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Fig. 11. W[ 12 CO]/W[ 13 CO] versus W[ 13 CO] 

The decrease of jy[ 12 CO]/W[ 13 CO] versus W[ 13 CO] 
(Fig. [HJ is caused by the saturation of the CO line. The few 
high ratio points at low W[ 13 CO] are attributed to the case 
in which both isotopes are optically thin and the asymptotic 
value of about 3 is reached when 13 CO becomes optically thick. 
These small ratios are found in the central parts of the cloud. 

The mean value of the ratio W[ 13 CO(l-0)]/W[C 18 O(l-0)] 
is 7.2±0.3. Within a central region of some 10' x 15' the ratio 
decreases, presumably due to 13 CO saturation, below the nor- 
mally assumed isotopic abundance ratio of 5.5. On the other 
hand, everywhere in the outer region of the area covered by 
C 18 observations the ratio increases to almost 10. Although 
13 CO emission is still not completely optically thin, the ob- 
served variation also indicates a change in the relative abun- 
dance of the two species. The 13 CO/C 18 ratio is enhanced to 
values >5.5 by 13 CO fractionation. The average ratio for the 
corresponding J =2-1 lines is 6. 8 ±2. 2. 

5. Column density distribution and mass of 
L1642 

5.1. LTE-analysis 

In order to calculate the column density, one must first de- 
termine the excitation temperature throughout the cloud. The 



Error bars reflect dispersion in each distance bin 



standard methods assume LTE conditions, and either 12 CO is 
assumed to be optically thick, or the optical depth ratio of two 
isotopic lines is assumed equal to the terrestrial abundance ratio 
(the assumed 12 CO to 13 CO abundance ratio is ~89). In addi- 
tion, an identical excitation temperature is usually assumed for 
both isotopes. 

In our case, both methods agree up to distances of ~30' 
from the cloud centre (the adopted cloud center is: -5', +10'), 
after which the 13 CO antenna temperature values are affected 
by noise. 

The excitation temperature has a nearly constant value from 
cloud centre to edge (Fig. fl3. 

In the central part of the cloud, 12 CO and 13 CO emission 
lines are certainly optically thick, and hence the derived tem- 
perature may not be representative of the bulk of the cloud. The 
best temperature estimate for the cloud centre is obtained from 
the 13 CO/C 18 ratio which gives a value of T ex =8 K. Based on 
the optically thick 12 CO emission, the mean temperature esti- 
mate isT ex =10.7K. 

The column densities and the cloud mass were estimated 
based on the 13 CO data. The estimates are derived from the 
J=l-0 line, since according to radiative transfer models (see 
Appendix lB.2l» this gives more reliable estimates than the J=2- 
1 transition. 

For LTE calculations, we assume a constant excitation tem- 
perature of 10K. Figure [Q] shows the H 2 column density map 
derived from 13 CO(1-0). The adopted conversion factor is 
N(H 2 )/ N( 13 CO)=l 10 6 . IRAS 04237-1419 is close to the col- 
umn density maximum, while IRAS 04325-1419 is located half 
way to the edge of the central core. The H2 column density was 
averaged in concentric rings of angular radius r from the centre, 
and the resulting radial column density distribution is shown 
in Fig. [14| The vertical error bars represent the dispersion of 
the values averaged at each radius. An r _1 distribution is ob- 
tained in the outer parts of the 13 CO cloud (r >5.6'), but the 
distribution in the core (r <5.6') is significantly flatter. Such 
behaviour is commonly observed in cloud cores (Bacmann et 
al. l2000l Ward-Thompson et al.. [T994l[T999t . and is in good 
agreement with both the previous 13 CO observations of LI 642 
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Fig. 13. Map of H2 column density derived from 13 CO(1-0) 
lines assuming LTE. The contours range from 1.0 10 20 to 5.6 
10 21 cm -2 by 5.0 10 20 cm -2 . The filled square and the triangle 
indicate the positions of IRAS 04237-1419 and IRAS 04325- 
1419 



ft 



Angular Radius i 



Fig. 14. The profile of the H2 column density derived from 
13 CO(1-0) observations. The dotted lines are for r _1 trends and 
dashed lines for r -0,2 . The centre is at Aa=-5', Ac^+10' (see 
Fig.O 

(Liliestrom. 119911 see her Fig. 7) and the density profile of dust 
(Laureiiset al. [T987b . 

The cloud mass can be calculated according to 



M(M ) 



5.74 10" 20 d\ mpc J J N(H 2 )cosSdSda, (4) 



where <iioop C =^/100pc, with d the distance of the cloud. Within 
the boundary T^( 13 CO(l - 0)) > IK, corresponding to 
if (H 2 ) > 6 x 10 20 cm -2 , we obtain from the iV(H 2 ) map as 
shown in Fig.lT3lM = (30 ± 10) M d 2 100pc - 59 M (assum- 
ing d=140 pc). 

5.2. Spherically symmetric non-LTE models 

In reality, the excitation temperature is not likely to be con- 
stant over the whole line of sight. This may lead to erroneous 
results in the normal LTE analysis, especially at low tempera- 
tures (e.g. Padoan et al. 2000). For this reason, we compared 
the LTE results, as described above, with non-LTE radiative 



transfer calculations. The comparison gives an estimate for the 
uncertainty of the results. Precise determination of the cloud 
mass requires detailed knowledge of the density and tempera- 
ture distributions. 

The 13 CO observations were modelled with a spherically 
symmetric cloud with radial density distribution n ~ r -2,0 , ex- 
cept for the central 10% of the cloud radius for which constant 
density was adopted. The density contrast between the centre 
and the cloud surface was set to 20. The radiative transfer prob- 
lem was solved with a Monte Carlo method (Juvela 1997}, and 
the results were compared with observed spectra which were 
averaged over concentric annuli. The density and the size of 
the model cloud were varied in order to find parameter values 
that best reproduced the observations. A value of 1 .0- 10 -6 was 
assumed for the fractional abundance of 13 CO. 

Isothermal models with T^ n between 8 and 14 K were first 
considered. The \ 2 values were practically independent of the 
temperature and it is clear that the temperature cannot be well 
constrained based on 13 CO observations alone. The mass es- 
timates ranged from 82 M at 8 K to -65 M at 14 K. The 
intensity ratio between the 2-1 and 1-0 lines was generally too 
low and, as expected for a microturbulent model, the line pro- 
files toward the centre positions were strongly self-absorbed. 

We further studied a set of models in which the kinetic tem- 
perature increase linearly with the radius. These did not lead 
to significantly better fits but the mass estimates were slightly 
lower and close to the LTE values. For a model where T^ n 
increases from 9K in the centre to 14 K on the cloud sur- 
face, the estimated mass was 61 Mq and the central density 
was 1.2- 10 4 cm 3 . For this model the fitted spectra are shown in 

Fig.na 

As can be seen the spherically symmetric models do not 
provide a good fit to the observations. At the centre position the 
observed spectra are clearly asymmetric and when observations 
are averaged over concentric annuli the presence of different 
velocity components or gradients leads to artificial line broad- 
ening or even multiple peaks in the averaged profiles. Such 
effects can be taken into account only in three-dimensional 
models (see Sect. 15.31 add Appendix iBl. The main shortcom- 
ing of the spherically symmetric models is, however, the low 
line ratio between the J=2-l and J=l-0 transitions. This could 
indicate either a much higher kinetic temperature, significant 
foreground absorption or more likely density inhomogeneities 
inside the source. 

The observed average intensity ratio 13 CO(2-l)/ 13 CO(l- 
0) is 0.69 (the average of observations with W( 13 CO(l- 

0) )>1 Kkms -1 ). Assuming LTE conditions, constant excita- 
tion along the line of sight, and an excitation temperature of 
T ex =10K (Tig.[T2l». such a ratio would imply an optical depth 
of 13 CO(2-l)^2. The 13 CO molecule is, however, not thermal- 
ized. In our calculations, the excitation temperature of the 2-1 
transition is always clearly below that of the transition 1-0. For 
example, in the centre of the model cloud with Tki n =12K the 
difference is almost three degrees. At higher densities 13 CO 
would become less subthermally excited and the ratio 13 CO(2- 

1) / 13 CO(1-0) would increase. Therefore, the observed high line 
ratios may indicate that the cloud is inhomogeneous and 13 CO 
is partly thermalized in spite of the low average density. 
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Fig. 15. Spectra from a spherically symmetric model cloud 
(solid lines) and observed spectra (histograms) averaged 
over different annuli centred at the position 4h35ml3.6s - 
14°15 / 28 // . In the model the kinetic temperature increased lin- 
early from 9 K in the centre to 14 K on the cloud surface. The 
radius of the annulus is given in each frame. Thick lines repre- 
sent 13 CO(1-0) and thin lines 13 CO(2-l) spectra. 



5.3. Three-dimensional non-LTE models 

The line profiles from a real, turbulent cloud are seldom as 
clearly self-absorbed as in the micro turbulent models. The 
cloud was therefore also modelled with fully three-dimensional 
models. In these isothermal models, the general density distri- 
bution was based on the LTE column density maps. Details of 
the calculations are given in Appendix lB.il The best fit was ob- 
tained at a significantly lower temperature than for the spheri- 
cal homogeneous models, Tki n =8 K, leading to a mass estimate 
of 70 Mq . However, the temperature cannot be accurately de- 
termined, and for the model with Tki n =14K, the average line 
ratio, 13 CO(2-1)/ 13 CO(1-0)=0.63, is still reasonably close to 
the observed value of 0.69. For any given temperature, the pre- 
dicted mass values were only slightly lower than in the spheri- 
cal models. 

We applied a standard LTE analysis to the spectra produced 
by the 3D models. This analysis showed that, in the parameter 
range considered here, the LTE column density estimates based 
on the 1-0 are reliable (see Appendix IB. 21 . The mass estimate 
from the LTE-calculations are indeed consistent with the 3-D 
modelling, when a kinetic temperature of T^ n ~ 12 K is as- 
sumed. 



where cr gas and <j v i r are the velocity dispersion of the gas de- 
rived from observation and the virialized velocity dispersion of 
the mean gas particle, respectively. The effective radius of the 
cloud is i? e ff = 0.97 dioo pc pc (R e f? = \foh\ a=semi-major 
axis, 6=semi-minor axis). A value of cr gas =0.17kms -1 is ob- 
tained from the 13 CO line width (corrected for instrumental 
and line opacity broadening). Adopting Tki n =12K and assum- 
ing <i=140pc, we find cr v i r) i_di m =0.21±0.02 and hence a ratio 
of ^^=0.81 ±0.08 for a homogeneous sphere. For a centrally 
condensed sphere with n oc r _1 , we obtain cr v i r =0.19±0.02 
and ^^=0.89±0.08. The values are very close to 1, suggest- 
ing that the cloud is probably in virial equilibrium. Wed note, 
that 13 CO traces only some 80% of the total cloud mass (e.g. 
Falgarone and Puget, 1988 ), while another 20% is in more dif- 
fuse cloud medium. 



6. Conclusions 

We have presented 12 CO, 13 CO and C 18 J=l-0 and J=2-l 
line observations of the translucent high latitude cloud LI 642. 
The main results of this study are summarized as follows: 

- Besides the conventional channel maps, we have also ap- 
plied the Principal Component Analysis (PCA) and Positive 
Matrix Factorization (PMF). The PMF method was found to 
produce an easily understandable, compressed presentation of 
the data that preserved all the important morphological and 
kinematic features of the cloud. 

- The morphological and kinematic analysis has shown a 
main structure at a velocity of 0.2kms _1 with an incomplete 
ring structure around it at larger velocities. Fainter emission 
extends to the north with velocities up to 1.6 km s _1 . 

This structure suggests an incomplete expanding shell, the 
front side of which is seen in the main component. 

An elongated morphology in the line of sight direction is 
also suggested by the kinematics. 

- In C 18 0, only two spatially well separated clumps can be 
distinguished. These structures do not coincide with the two 
embedded IRAS sources, which are offset by 5-7' from the 
main peak. 

- We confirm that LI 642 follows an r _1 density distribution 
in its outer parts and a flatter distribution with n oc r -0,2 for 
the core. The cloud is probably in virial equilibrium. 

- The LTE estimate of the cloud mass is 59 Mq. 

- We performed three-dimensional radiative transfer mod- 
elling of the 13 CO observations and obtained a cloud mass 
estimate of ~10Mq. The cloud temperature is, however, not 
well constrained, leading to large uncertainties in the mass es- 
timates. 



5.4. Dynamical state of the cloud core 

We investigate the dynamical state of LI 642, following the 
method by Liljestrom (I1991L who estimated the ratio 
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40 20 -2G-+0 

Fig. A.l. 12 CO(1-0) channel maps. The velocity (in km s -1 ) is indicated on each panel. The contour interval is 0.3 K km s -1 and 
the first contour is at 0.2 K km s -1 

Appendix A: Alternative methods for the study 
of cloud kinematics 

A.1. Channel maps 

The 12 CO, 13 CO and C 18 channel maps are shown in Figs. 
IA.1I to IA.4I In 12 CO, the first significant emission appears 
at -0.5 km s -1 around position (-3', +10') and, as the veloc- 
ity increases, emission extends towards both south and west. 
Above -0.1 kms -1 , the general morphology remains the same 
with only small additional patches appearing. These suggest 
the presence of small internal motions and clumpy structure. 
Finally, above 0.5 kms -1 , the cloud breaks up into two sepa- 
rate structures. The first one (structure A; see Fig. IA.2l) is elon- 
gated in the east- west direction and is located around Ay ~20'. 
The second one (structure B) extends from the centre towards 
the south. 

Structure A remains smooth between 0.6 and 0.9 kms -1 . 
At higher velocity, it fragments with persistent clumps at (-20', 
+200 and (+32', +15'). The latter is well identified in 12 CO(l- 
0) but is less clear in 12 CO(2-l). The velocity of these clumps 
is ~1.2kms _1 . Structure B fades after 1 kms -1 , the mean ve- 
locity being ~0.8 km s -1 . At higher velocities, the main struc- 
ture is a faint elongated emission emanating from the (-20', 
+200 clump. 

Fig. IA.5I shows the velocity diagram for 12 CO(1-0) ob- 
servations along the line between positions (40', -400 and (- 
40' ,400- The average velocity increases from ^0.3 kms -1 in 
the south to ~1 kms -1 at the northern end where the velocity 
increases all the way to the edge of the map. The three preem- 
inent clumps at ~(0', 00, (-10', +100 and (-20', +200 are vis- 
ible as separate maxima. However, only the last one (at offset 




Velocity (km/s) 



Fig.A.5. Velocity diagram of 12 CO(1-0) observations along 
the line from position (40 / ,-400 to (-40',400. 



30' in Fig. IA.5l) is clearly separated from the main cloud based 
on its radial velocity. The gradual change between offsets 15' 
and 30' does not necessarily require a continuous velocity gra- 
dient but could be caused by the superposition of emission at 
two separate velocities. The range of velocities is large, both at 
the location of the northern clump at (-20', +200 (offset +300 
and in the south (offset -150- However, the largest dispersion 
is found at offset 7 f where the emission extends to high posi- 
tive radial velocities, especially compared with the velocity of 
the main emission. The location is some 5' from IRAS 04325- 
1419, the source of a known outflow (Liljestrom et al. 1989). 
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40 20 -2(^40 

Fig. A.2. 12 CO(2-l) channel maps. The velocity (in km s _1 ) is indicated on each panel. The contour interval is 0.2 K km s _1 and 
first contour is at 0.2 K km s _1 . On this figure, structures A and B, mentioned in the text, are indicated on the 0.6 km s _1 channel 



The 13 CO emission is less extended, but it exhibits simi- 
lar kinematics as 12 CO. In 13 CO(2-l), map emission clumps 
are clearly seen around (+6', +120, (-6', +11.50 and (+9', -20- 
The (-20', +200 emission clump is still dominant, with a mean 
velocity at around 0.8 km s _1 . 

C 18 probes the deepest and densest parts of LI 642. 
Compared with the 12 CO and 13 CO lines, it is spatially much 
more concentrated and is restricted to a velocity range be- 
tween -0.3 and 0.5 km s _1 . The emission mostly occurs around 
— 0.15kms _1 . Two main parts can be distinguished around 
Ax —5' (part A) and Ax ~3' (part B) (see Fig. [Oil . Both 



have several patches appearing at different velocities: in A at (- 
4', +90 04r —0.1 km s- 1 ), (-6', +90 (Msr -0.1 km s" 1 ), (-4', 
+60 04r -0.15 km s- 1 ) and (-4', +3.50 04r -0.25 km s" 1 ), 
and in part B at (+2', +60 (Msr -0.1 kms" 1 ) and (+6',+90 
(Vi S r —0.3 km s _1 ). These could either be independent clumps 
or velocity cuts of a continuous density distribution where there 
are gas motions along the line of sight. The latter hypothesis is 
more likely, because the features do not have a clear counter- 
part on peak intensity and line area maps. 
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Fig. A.3. 13 CO(1-0) (upper) and 13 CO(2-l) (lower) channel maps. The velocity (in kms -1 ) is indicated on each panel. The 
contour interval is 0.1 Kkm s _1 for 1-0 and 0.05 Kkm s _1 for 2-1, both starting at 0.1 Kkm s _1 



A.2. Principal component analysis 

The spectral data consist of an ensemble of n objects (profiles) 
each with p attributes (intensity at each velocity channel) and 
can thus be formally represented as an n x p matrix X. The 
projection of the data onto an arbitrary axis u a is given by y a 
= Xu a . The principal component analysis (PCA; Murtagh & 
Heck 1 19871 determines a set of orthogonal axes u a , such that 



the projection of X upon each consecutive u a always max- 
imises the variance along that axis. Mathematically, the base 
is determined by computing the eigenvectors u a of the covari- 
ance matrix. The method gives an objective way of determining 
the characteristic radial velocities of structures which may not 
even be discernible upon visual inspection. 
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Fig. A.4. C 18 O(l-0) (upper) and C 18 0(2-l) (lower) channel maps. The velocity (in kms x ) is indicated on each panel. The 
contour interval is 0.03 K kms -1 for 1-0 and 0.02 K kms -1 for 2-1, the lowest contour lying at 0.1 Kkms -1 



The relative importance, C a , of the axis a is defined as 
C a — v^ Aq \ , i.e. the eigenvalue X a indicates the percentage 

of the variance of the observed variables explained by the cor- 
responding principal component. Usually two or three factorial 
axes with the highest C a value are enough to describe the data. 
Further details about the application of PCA to the study of 
the interstellar medium can be found in Heyer and Schloerb 
( 1997 ), Abraham et al. <l2QQ0b and Ungerechts et al. <fT997t . 

PCA was applied independently to 12 CO(1-0), 12 CO(2-l), 
13 CO(1-0), 13 CO(2-l), C 18 O(l-0) and C 18 0(2-l) data cubes. 
For 13 CO and C 18 lines, the eigenvectors are very irregular. 
Normal PCA analysis does not take into account error esti- 
mates, and a good signal-to-noise ratio is important, especially 
as PCA is computed using normalised profiles. Therefore, we 
limited the PCA analysis to the 12 CO(1-0) data. 

From 12 CO(1-0) data a 648x40 matrix was constructed 
containing 648 observed spectra with 40 velocity channels. The 
velocity range is [-1.6,+2.6]kms _1 with a channel width of 
0.1 km s _1 . Before PCA analysis, the data are standardised, i.e. 
each spectrum is mean-centred and normalised with its stan- 
dard deviation. 

Figures lA!6l and Ia771 show our results. The principal com- 
ponents 1-4 account for 71%, 13%, 7% and 3% of the total 
variance, respectively. Together, these represent 95% of the to- 
tal variance in our data. The remaining components are mostly 
pure noise and are discarded. The importance of the first com- 
ponent shows that the spectra are not independent from each 
other. As expected, there is a rather continuous variation of line 
morphology throughout the cloud. 



The first principal component, PCI, represents the aver- 
age spectral profile and the associated map shows the clumpy 
distribution of the main cloud. The emission corresponds to 
the structures seen in channel maps at ~0.5kms -1 . The first 
eigenvector is roughly gaussian, with position ~0.6 km s _1 and 
width 1.4 km s _1 . The following PCs describe successive cor- 
rections to the first approximation given by PCI . Looking at the 
profiles, PC2 suggests a shift in velocity and PC4 a change in 
the linewidth. PC3 indicates either a small change in the radial 
velocity or, more likely, the presence of some narrower spectra. 

PC2 mainly indicates velocity variations inside central re- 
gion of LI 642 and separates this from the northern part of the 
field. The eigenvector indicates that these sub- structures are 
around 0.2 km s _1 while the northern part has a mean velocity 
of 1.2 kms -1 . PC3 indicates a region of emission around (0', 
00 and (-9', +8.50 with a velocity around 0.8 kms -1 . Finally, 
PC4 shows substructures in the core at velocity ^O.lkrns -1 
especially around (-2.5', +90, (-8, +170 and (-8.5', +90- The 
map of PC4 resembles the channel maps between 1.5 and 
1.8 km s -1 . The northern part appears patchy, but it is difficult 
to say if the structure is caused by real clumps or is simply due 
to noise fluctuations. 

The PCA analysis follows a hierarchical scheme: first, it 
shows the central part of LI 642 (0.6 kms -1 ), secondly, it dis- 
criminates this from the northern part (1.2 kms -1 ). Next it 
reveals additional velocity structure inside LI 642 itself (0.2, 
0.8 kms -1 ), and finally fainter emission at more positive ve- 
locity in the extreme north (~1.5 km s -1 ). 
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Fig. A.6. Eigenimages of the four first principal components Fig. A.7. Eigenvectors used to generate the eigenimages shown 
calculated from the 12 CO(1-0) data. Negative contours are in Fig. lA^I 
dashed (contour spacing 0.1), grey scale represents positive val- 
ues (for PCI contour spacing is 0.05 else it is 0.2) and the full 
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A3. Comparison of the analysis methods 

Channel maps are the standard method for analysing kinematic 
information of molecular line maps. The presentation is mostly 
limited by the noise. If the channel width is kept constant, the 
noise will dominate in the line wings. If, on the other hand, 
noise is suppressed by integration over larger velocity inter- 
vals, one loses all information of smaller velocity structures. 
Furthermore, the signal level varies significantly over the map, 
and the selection of integration intervals and contour levels is 
necessarily a compromise. However, one can visually recog- 
nise continuous regions of emission even when noise is quite 
significant. Channel maps do not provide direct information on 
the spectral profiles because each channel is treated completely 
separately. This information is obtained indirectly by looking 
at the maps of consecutive velocity channels. 

Both PCA and PMF present the observations as a linear 
combination of basic spectral components which are extracted 
from the data. The methods carry out a global fit, and can there- 
fore discern features that are too weak to be detected in any in- 
dividual spectrum but which are repeated in many observations. 
When the intensities of the spectral components are determined 
for a given spectrum, the fit is done using all the relevant chan- 
nels. The methods are therefore most useful in the case of spec- 
tra with low signal-to-noise ratios. The spectral features cannot, 
of course, be separated unless their relative intensities change 
over the map. PCA and PMF decompositions compress the in- 
formation contained in the observation matrix to be presented 
with just a few components. The result should therefore be eas- 
ier to interpret. In the case of PCA, however, this is complicated 
by negative values present both in the profiles and in the com- 
ponent maps. 

In the case of our LI 642 observations, the three methods 
agree in the main features. Comparison of the different meth- 
ods shows that the PMF method is especially able to produce 
a presentation that is at the same time compact and easy to 
interpret. For example, a PMF decomposition with r=4 repro- 
duces practically all of the structure visible in the very complex 
channel maps of 12 CO. Because of the lower signal-to-noise ra- 
tio, only three components could be detected in the C 18 data. 
When the rank of the factorization is lower than the rank of the 
observation matrix, the presentation is not exact. However, if 
r is sufficiently large, the fit residuals mostly consist of noise. 
The study of residuals of the PCA or PMF fits is a very power- 
ful tool for finding deviating spectra (see Sect l3.ll) . Residuals 
identify abnormal spectra and show how their profiles differ 
from the average behaviour. 

In channel maps, one has the tendency to interpret inten- 
sity variations of separate components as continuous velocity 
gradients. PMF and PCA explicitly assume that the emission 
is the sum of components at fixed radial velocities. Velocity 
gradients are approximated as the sum of two or more spectral 
components at different velocities, and with gradual changes in 
the weights of the components. In the northern part of LI 642, 
the channel maps of 12 CO J=l-0 and J=2-l suggest a small 
velocity gradient towards the north. PMF decomposition with 
r=4 presents the northern part with just one velocity compo- 
nent. With r=5, the emission is divided into two components, 



suggesting the presence of a velocity gradient. With higher r, 
the emission is not divided any further. Within the noise limits 
the observations are, therefore, consistent with both a small ve- 
locity gradient and the presence of two separate emission com- 
ponents. 

Appendix B: Three-dimensional radiative 
transfer model 

B.1. Mass estimates 

The 13 CO J=l-0 and J=2-l observations of L1642 were mod- 
elled with a fully three-dimensional model that consisted of 46 3 
cells arranged on a Cartesian grid. Along two dimensions, the 
column density distribution follows the distribution obtained 
from the LTE calculations, and this fixes the size of the model. 
Along the third dimension (i.e. along the line-of- sight), the den- 
sity distribution is assumed to be gaussian, with the extent of 
the cloud similar along the line-of- sight as in the plane of the 
sky. In the model, the linewidths are caused in roughly equal 
parts by microturbulence within cells and macroscopic turbu- 
lence between cells. The mean velocity of cells on each line- 
of-sight was shifted according to the velocity of the observed 
13 CO(1-0) line. This enables direct comparison with the ob- 
served spectra, but also creates velocity gradients that have 
a small effect on the excitation. The density distribution was 
modified according to a volume filling factor, /, by reducing 
the density by two orders of magnitude in a randomly selected 
fraction l-f of the cells. This is only a crude approximation 
of the true cloud structure, but such a model has already been 
shown to reproduce some of the effects caused by density inho- 
mogeneities (see e.g. Park, Hong, & Minh 1996] Juvela|l997). 
The three-dimensional model is fundamentally different from 
the one-dimensional models discussed in Sect. 13 and the com- 
parison of the results should give a good idea of the overall 
uncertainty of the mass estimates. 

The parameter / was varied between 0.1 and 1.0, and 
model calculations were carried out for Tki n =8-14K. The best 
fit with observations was obtained with Tki n =8 K, / =0.4 and 
peak density of 1.7xl0 4 cm -3 . The total mass of the model 
cloud was 70 Mq, very similar to the values obtained with 
spherically symmetric models. The x 2 -value of the fit was 
less than 10% higher at 12 K suggesting that the temperature 
is not well constrained. Since the cloud size was fixed, the 
mass decreases when the assumed temperature increases: for 
^kin=14K it is 37 Mq. The volume filling factor of the best 
fitting models increases with temperature, and was already 1.0 
at 14 K. The /-parameter had surprisingly little effect on the 
results, particularly at higher temperatures. However, because 
of the macroturbulence, all models contained equal amount of 
inhomogeneity in the velocity space. 

B.2. Comparison of true and LTE column densities 

In order to further investigate the errors introduced by the LTE 
assumption, we computed LTE column density estimates based 
on 3D cloud models with Tki n =10K and volume filling fac- 
tors /=1.0 and /=0.1. The model parameters were adjusted 
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so that the model approximately reproduced the 13 CO obser- 
vations of LI 642. Spectra were calculated for the 12 CO(1-0) 
and 13 CO(1-0) lines, and LTE analysis of the line ratios of 
the calculated lines resulted in average excitation temperatures 
of T ex =8.9 K for the model with f=l.O and T ex =8.4 K for 
the model with /=0.1. LTE column densities were computed 
based on these average excitation temperatures and the calcu- 
lated 13 CO(1-0) spectra. In Fig. fly] the LTE column densi- 
ties are plotted against the known, true column densities in the 
models. 

For column densities derived from the 13 CO(1-0) line, the 
LTE calculations are quite accurate. In the case of f=0A, the 
scatter is larger because of the inhomogeneity of the cloud, and 
at high column densities the LTE-estimates may be slightly bi- 
ased towards low values. The LTE estimate is not very sensi- 
tive to the assumed excitation temperature, and even a change 
of 2 K in the value of T ex has little effect on the derived val- 
ues. The volume-averaged, true excitation temperatures ob- 
tained from the radiative transfer models are close to 5 K, i.e. 
clearly below the LTE estimates of 8-9 K. However, the large 
difference is mostly caused by low-density gas at the edges of 
the cloud, and the bulk of emission comes from denser regions 
which are much closer to thermalization. 

We also plot in Fig. lB.ll LTE estimates of the column den- 
sity derived from the model-calculated J =2-1 lines, assuming 
excitation temperatures that are 1 K lower than for the J =1- 
line. The column density values are seen to be too low by 
up to a factor of two. In the case of /=0.1, the density of the 
remaining clumps is higher, and since the gas is closer to ther- 
malization, the LTE estimates are also closer to the true column 
density values. The results are more sensitive to the value of 
T ex than in the case of the J =1-0 line. The column density 
values increase if one assumes a lower excitation temperature. 
However, the LTE estimates of the excitation temperature are 
less than 1 K lower than for the J =1-0 transition, and the 
same applies to the average excitation temperature in the nu- 
merical models. The results show that the column density es- 
timates based on the J =2-1 transition are inherently more 
uncertain than those derived from the J =1-0 observations. 
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